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Abstract 

We study stimulated Raman scattering (SRS) in laser-fusion conditions with the 
Eulerian Vlasov code ELVIS. Back SRS from homogeneous plasmas occurs in sub- 
picosecond bursts and far exceeds linear theory. Forward SRS and re-scatter of back 
SRS are also observed. The plasma wave frequency downshifts from the linear dis- 
persion curve, and the electron distribution shows flattening. This is consistent with 
trapping and reduces the Landau damping. There is some acoustic (cu oc k) activity 
and possibly electron acoustic scatter. Kinetic ions do not affect SRS for early times 
but suppress it later on. SRS from inhomogeneous plasmas exhibits a kinetic enhance- 
ment for long density scale lengths. More scattering results when the pump propagates 
to higher as opposed to lower density. 

1 Introduction and Code Model 

Laser-plasma interactions must be controlled for inertial fusion to succeed. This paper 
examines stimulated Raman scattering (SRS), or the parametric coupling of a pump light 
wave (the laser, mode 0) to a daughter light wave (mode 1) and an electron plasma wave 
(EPW, mode 2). Kinetic effects, such as electron trapping, in the daughter EPW are seen 
to be important in back SRS (BSRS). 1-D kinetic simulations presented here show BSRS 
much greater than coupled- mode theory for both homogeneous and inhomogeneous plasmas. 
Strong nonlinearity and non-thermal electron distributions are seen to result. 

ELVIS 1 is a 1-D Eulerian Vlasov code that evolves the distribution function f s (s = 
species; e for electrons) on a fixed phase-space grid. It uses operator splitting for the time 
advance [2], HH and cubic spline interpolation for shifting f s in position (x) and momentum 
(p x )- Light waves are linearly polarized in y. The ions can be immobile or kinetic. The 
governing equations are 

[dt + (p x /m s )d x + (Z s e) (E x + v ys B z ) d p J f s = v Ks (x) (n s f 0s - (1) 
d x E x — — Z s n s m s d t v ys = eZ s E y (2) 
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Figure 1: (a) Reflectivity for homogeneous run with immobile (solid curve) and kinetic ions 
(dotted curve) . (b) Spectrum of reflected light for immobile ions. "Fre" and "SEAS?" label 
BSRS re-scatter and possible electron acoustic scattering. 

(d t ±cd x )E ± = -~Y Z s n s v ys E ± =E y ±cB z (3) 

e ° ; 

A number-conserving Krook relaxation operator is included, with relaxation rate vk s (x) and 
equilibrium Maxwellian /o s (/ dp fo s = 1). We use a large vks ~ 0.2u> p {uit = nae 2 / (eom e )) 
at the edges of the finite density profile to absorb plasma waves generated by SRS and 
prevent reflection. A nonzero central value of vks can also mimic sideloss from a laser 
speckle. We advance E ± without dispersion by shifts of one x gridpoint, which imposes 
dx = c dt. 

2 Simulations Results 

We simulate a pump laser with Ao = 351 nm (vacuum) and intensity Iq = 2 x 10 15 W/cm 2 
impinging from the left (E + contains the pump) on a finite plasma with a flat central region 
75.1 nm wide of density no = 0.1n c (critical density n c — nouiy/uip) and temperature T e = 3 
keV. Since Vlasov codes are low-noise there are no numerical fluctuations for SRS to grow 
from. We therefore inject a counter-propagating seed light wave via E~ with Ai a = 574 nm 
and I\ = 10 _5 /o- This light has the maximum linear BSRS growth rate and couples to an 
EPW with fc2A_o = 0.357 and a Landau damping rate V2 — 0.038w p (Xd = VTe/^ p , = 
T e /m e ). The x grid spacing is dx = 11.9 nm, our algorithm requires dt = dx/c, and we use 
a p x grid spacing dp = 0.0437wTe"ie- SRS is convective for these parameters with a spatial 
gain rate a = 0.019 /im -1 , giving a linear reflectivity i?jj n = 1.72 x 10~ 4 . The numerical R, 
shown as the solid curve in Fig. 1(a), is well above this level. R comes in sub-picosecond 
bursts and has a time average from 1 ps to the run end of R av = 13.8%. vks ^ only at the 
edges. Repeating the run with a nonzero central vks shows a sharp cutoff of the reflectivity 
for vks ^ 10~ 3 w p . 

The streaked spectrum of reflected light E~ at the left edge is shown in Fig. 1(b). 
Almost all the energy is contained in BSRS. u>i s — 1.93o-> p is the seed light frequency. 
Initially BSRS occurs at uj\ s but upshifts for t > 2 ps, corresponding to a downshift in u>2- 
The weak signal near lu p labeled "Fre" is the forward Raman re-scatter of upshifted BSRS 
light. The longitudinal E x spectrum in Fig. 2(a) reveals the plasma wave at the matching k 
and lo. Re-scatter is only possible due to the upshift in u>\, since u p > w\ s /2. The features 
slightly below luq and slightly above 2.5uj p , labeled "SEAS?" , may be related to scattering off 
the acoustic longitudinal activity discussed below 0]. The transmitted light (E + ) spectrum 
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Figure 2: (a) E x (k,uj) spectrum for homogeneous, immobile ions run. "Fre", "FSRS", and 
"ac" label BSRS re-scatter, FSRS, and acoustic activity. The black dots are the linear EPW 
dispersion curve, (b) rms-averaged E x (x,£). 

(not presented here) shows for t >3 ps weak levels of both FSRS and the anti-Stokes line 
(u> = luq + u>p) of the pump, even though neither of these is externally seeded. 

The longitudinal electric field spectrum E x (k,u>), depicted in Fig. 2(a), reveals that the 
BSRS plasma-wave activity is downshifted in frequency from the linear dispersion curve. 
This is qualitatively consistent with the frequency downshift due to electron trapping 0, 
and periods of larger downshift correspond to larger EPW amplitude. The downshifted 
EPW connects with an acoustic- like feature (u) oc k) that extends to uj = 0. There is 
another, lower phase velocity acoustic streak, strongest for uj < 0.2w p . Weak scattering 
off them may account for the "SEAS?" features in Fig. 1(b). In addition, plasma waves 
corresponding to FSRS and re-scatter of BSRS occur on the EPW dispersion curve. Fig. 
2(b) presents the t and x rms-averaged E x (x,t), which shows the EPWs occur as a series 
of wide pulses that move parallel to the laser (that is, to the right). The group velocity 
matches the slope of the BSRS plasma waves. Near the left edge some pulses propagate 
opposite the laser. 

The electron distribution f e forms phase-space vortices at the EPW phase velocity v P 2 = 
(=0.264c for the linear EPW). The space-averaged (/ e ), displayed in Fig. 3(a), is 
flattened due to trapping in this region. Landau damping (~ df/dv) is greatly reduced by 
flattening, which may thereby enhance the reflectivity 6 . When the EPW amplitude is 
large f e is quite flat, and only for brief periods (< 0.1 ps) do we see a small bump (region 
of d{f e )/dp x > 0) form slightly above v P 2- 

The run was repeated with kinetic helium ions (m, = 4m p , 7j = 750 eV) and yielded 
the dotted R in Fig. 1(a). Early in time R is the same for immobile and kinetic ions, while 
later in time they diverge. For the last 2 ps the reflectivity is very low with kinetic ions. 
We do not see evidence of Langmuir decay instability (EPW — ► EPW + IAW (ion acoustic 
wave)) of the EPW or stimulated Brillouin scattering (photon — > photon + IAW) of the 
pump. Instead, very high E x activity develops on the left edge of the box around t = 5 ps, 
involving large ion density fluctuations; BSRS is minimal after this. Most of the spectral 
power in this activity is concentrated in the FSRS and re-scatter of BSRS features. Further 
study of the role of ions is underway. 

In an inhomogeneous medium, the k matching condition fco = k\ + &2 for a three-wave 
interaction can only be satisfied at one point. Away from this point the detuning limits the 
interaction. We performed ELVIS simulations for the same parameters as the homogeneous 
run with kinetic ions discussed above. However, the central region of the density profile now 
has a linear gradient extending for 100 fixa. We vary the endpoint densities and thereby 
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Figure 3: (a) Space-averaged f e over central 5.62 /im of box for homogeneous, immobile 
ions run. (b) Average reflectivity for inhomogeneous runs with L n = (167, 250, 500) /im. 
(Solid, dashed) curves are the coupled-mode (strong damping limit, Rosenbluth undamped) 
steady-state R. 

change the density scale length L n = n/{dn/dx). The reflectivities for several L n are 
shown in Fig. 3(b), for the pump propagating toward higher and lower densities. Also 
plotted is the steady-state R predicted from the coupled- mode equations, solved in the 
limit of strong damping for the EPW (solid curve) as well as the Rosenbluth undamped 
result (dashed curve). Both coupled- mode calculations give the same R for both directions 
of pump propagation, yet the simulation consistently shows higher R for /c ||Vn. We are 
formulating a theory to explain the high R and the role of pump propagation direction. 

3 Conclusions 

Vlasov simulations of SRS show strong enhancement of the scattering over coupled-mode 
predictions for both homogeneous and inhomogeneous plasmas. The resulting plasma waves 
do not satisfy the linear dispersion relation. The electron distribution shows large trapping 
and flattening. The role of sideloss and ions need to be further examined, and analytic 
models that explain these findings need to be developed. 

This work was supported by the US Dept. of Energy under grant No. DE-FG02- 
91ER54109 at MIT and under contract No. W-7405-ENG-48 at LLNL. 
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